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We calculate the properties of the two-band Hubbard model using the Dynamical Cluster Ap- 
proximation. The phase diagram resembles the generic phase diagram of the cuprates, showing a 
strong asymmetry with respect to electron and hole doped regimes, in agreement with experiment. 
Asymmetric features are also seen in one-particle spectral functions and in the charge, spin and d- 
wave pairing susceptibility functions. We address the possible reduction of the two-band model to a 
low-energy single-band one, as it was suggested by Zhang and Rice. Comparing the two-band Hub- 
bard model properties with the single-band Hubbard model ones, we have found similar low-energy 
physics provided that the next-nearest-neighbor hopping term t' has a significant vahie(t' /t m 0.3). 
The parameter t is the main culprit for the electron-hole asymmetry. However, a significant value 
of t' cannot be provided in a strict Zhang and Rice picture where the extra holes added into the 
system bind to the existing Cu holes forming local singlets. We notice that by considering approxi- 
mate singlet states, such as plaquette ones, reasonable values of t', which capture qualitatively the 
physics of the two-band model can be obtained. We conclude that a single-band t-t'-U Hubbard 
model captures the basic physics of the cuprates concerning superconductivity, antiferromagnetism, 
pseudogap and electron-hole asymmetry, but is not suitable for a quantitative analysis or to describe 
physical properties involving energy scales larger than about 0.5 eV. 



I. INTRODUCTION 

The theory of the cuprate high temperature supercon- 
ductors remains one of the most important and daunt- 
ing problems in condensed matter physics. The high 
T c cuprate superconductors are layered materials with 
relatively complex structures and chemical composition. 
They are highly correlated, with an effective bandwidth 
roughly equal to the effective local Coulomb interac- 
tion. The short-range correlations are known to play a 
paramount role in these materials. Therefore, the Dy- 
namical Cluster Approximation (DCA)i, which treats 
short-range correlations explicitly and the long-range 
physics at the mean-field level, is an ideal tool for the 
investigation of these systems. 

A common characteristic all the cuprate materials 
share is the presence of quasi-two-dimensional CuOi 
planes. These planes are commonly believed to contain 
the low-energy physics. However, the full complexity of 
the orbital chemistry of just the CuOi planes and the 
strong Coulomb repulsion on the Cu ions would lead to 
models which are very difficult to study with conventional 
techniques. 

The cuprates are characterized by a very rich, but also 
in many respects very intriguing physics. The undoped 
materials are antiferromagnetic (AFM) insulators with a 
gap of approximatively 2 eV. Upon doping the AFM is 
destroyed and the system becomes superconducting (SC). 
At small doping, in the proximity of the AFM phase, 



the normal state physics cannot be described in terms of 
Fermi liquid theory and is characterized by the presence 
of a pseudogap. An essential demand of every successful 
theory is to capture all these fundamental features at the 
same time. 

Experimental data show that the phase diagram and 
other physical characteristics like the density of states 
(DOS) near the Fermi level of the hole-doped and 
electron-doped materials are very differenliSiiii. There 
could be many reasons for this asymmetry. The electron 
and the hole doped materials are physically different and 
apart from the CuOi planes, they contain different el- 
ements and chemical structures. These structural and 
compositional differences can influence the low-energy 
physics. Therefore in this paper, we use DCA to ad- 
dress whether the physics of a pure CuOi plane contains 
this asymmetry or if the origin of the asymmetry in real 
materials comes from other influences. 

Different models for describing the physics of a CuOi 
plane were proposed by various authors. Photoemission 
experiments in the insulating parent material show that 
the first electron-removal states have primarily oxygen 
character; whereas, the first electron-addition states have 
d character, already suggesting a strong asymmetry. This 
places these materials in the charge-transfer gap region 
of the Zaanen-Sawatzky- Allen scheme 5 . Early on, con- 
sidering the ligand field symmetry and band structure 
calculations^*^, it was realized that the most important 
degrees of freedom are the Cu d x 2_ y 2 which couple with 
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the in-plane O p orbitals. Therefore, one of the first mod- 
els proposed to describe the physics of high T c materials 
was the so called three-band Hubbard model introduced 
by Varma et al& and Emery et ali^, which considers 
explicitly both the oxygen p a and the cooper d x i_ y i or- 
bitals. In fact, because the direct oxygen-oxygen hopping 
is neglected, only the combination of oxygen orbitals with 
x 2 — y 2 symmetry couples with the d orbitals, and the 
above proposed three-band model reduces to a two-band 
one. 

However Zhang and Rice (ZR) 1 ^ argue that the low- 
energy physics of the hole doped superconductors can be 
described by a single-band model. Starting from the two- 
band model, Zhang and Rice claim that an extra hole 
added into the oxygen band binds strongly with a hole 
on the Cu, forming an on-site singlet. This singlet state, 
which has zero spin can be thought as moving through 
the lattice like a hole in an antifcrromagnetic background. 
Consequently, the physics can be described by a one-band 
t-J model. 

Pertinent criticism to these simplified models were 
raised by various authors. With respect to Cu degrees 
of freedom, H. Eskes et alw^ stressed the possible impor- 
tance of the other d orbitals, showing that they should 
be explicitly considered when physics which implies exci- 
tations with energy larger than w 1 eV is involved. How- 
ever, these criticisms do not concern us for the present 
study since we are interested only in physics at energies 
lower than rj 0.5 eV . 

Investigating the relative importance of various param- 
eters describing the CuOi planes it was realized early 
on that in addition to the Cu on-site Coulomb repul- 
sion (Udd ~ 8 eV) and Cu-0 hopping integral (t p d ~ 
1.3 — 1.5 eV) the 0-0 hopping integrals result in a large 
O 2p band width (W s=s 5 eV), indicating that these 
should be included explicitly in any theoryi2iiii4. There- 
fore, using the DCA technique as a means of including all 
these most important parameters and bands, we address 
two major problems in this paper: The physics of the 
CuOi plane (including a detailed study of the electron- 
hole asymmetry) and the reduction of the multi-band 
model to a single-band model. 

Regarding the reduction to a one-band model, one of 
the most serious criticisms to ZR theory is the neglect of 
the O 2p band structur e - 15 ' 16 . The natural tendency of 
the finite oxygen bandwidth is to delocalize and desta- 
bilize the ZR singlets. The question arises whether the 
low-energy states (i.e. the ZR singlets) are still well sepa- 
rated from the higher energy states (i.e. the non-bonding 
oxygen states). Otherwise, the reduction to a single-band 
model which neglects these high energy states is not pos- 
sible. This problem was previously considered by Eskes 
and Sawatzkjii& within an impurity calculation approach, 
but there, unlike in the DCA approach, both the spatial 
correlation effects and the dispersion of the low-energy 
states were neglected. 

Another important objection to ZR theory was first 
raised by Emery and Reiter— , and regards the nature 



of the low-energy states. Are these states real singlets 
which can be mapped onto holes, or does the hole on 
the O bind into a more complicated state which involves 
more than one Cu hole? Choosing a particular solvable 
example, which considers the Cu spins arranged ferro- 
magnetically, they showed that the low-energy states are 
in fact an admixture of the Zhang-Rice singlets and the 
corresponding triplets. This implies a nonzero value for 
the oxygen spin, and destroys the equivalence of these 
states to holes. However, it is not clear if the situation is 
similar in the cuprates, i.e. if the ZR singlet-ZR triplet 
admixture is significant. But the merit of Emery and Rc- 
iter is to emphasize that the fact that, as a consequence 
of the strong Cu-0 hybridization low-energy states well 
separated from the non-bonding oxygen band states ap- 
pear, does not necessarily mean that the physics can be 
reduced to a single-band model. 

The third problem we address regarding the reduction 
to a single-band model is the estimation of the single- 
band parameters. We notice that different approxima- 
tions result in different values of the parameters. Espe- 
cially the magnitude of the next-nearest-neighbor hop- 
ping is very dependent of the initial assumptions. For 
example, if we assume that the hole addition low-energy 
states are genuine ZR singlets, i.e bound states between 
a Cu hole and a orthogonal Wannier oxygen orbital, we 
obtain a negligible next-nearest-neighbor hopping 1 —. On 
the other hand, if we consider the low-energy states to be 
plaquette singlets, i.e. bound states between a Cu hole 
and a hole on the state formed by the four oxygen around 
the Cu, the value of the next-nearest-neighbor hopping 
is significant^. Of course, because of the nonorthogonal- 
ity of the plaquette states, the plaquette singlets are not 
genuine singlets and therefore they cannot be rigorously 
mapped into holes. However, because their overlap with 
the local singlets is large (96%) 11 ' 17 , it is still possible 
that this approximation is good. 

Our calculations show that a multi-band model and 
a single-band t-t'-U Hubbard model with a significant 
value of the next-nearest-neighbor hopping exhibit a sim- 
ilar low-energy physics. The essential parameter needed 
for the agreement is the next-nearest-neighbor hopping, 
t' . This parameter is also the main culprit for the ob- 
served electron-hole asymmetry. However, as mentioned 
above, the large value of t' cannot be obtained in a strict 
ZR picture. Thus our results also implicitly indicate that 
the multi-band model cannot be rigorously reduced to a 
single-band model. Therefore, besides showing the simi- 
larities between the two models, we also point out their 
significant differences in this paper. 

The final conclusion is that a single-band t-t'-U Hub- 
bard model, with a significant value of t' , captures the 
basic physics of the cuprates, and thus is suitable to 
describe the AFM, pseudogap and SC physics together 
with the relevant asymmetries observed in the phase dia- 
gram, in the one-particle spectra and in the two-particle 
response functions. However, we believe that it is not 
suitable for a quantitative material specific analysis, for 
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describing the higher energy spectroscopic features as in 
optical spectroscopy or resonant inelastic x-ray scatter- 
ing, or for studying more subtle features related to the 
finite value of the spin on the oxygen. 

This paper is organized as follows. In Section [H] the 
two-band Hubbard model and the DCA technique is in- 
troduced. Our two-band model takes fully into account 
the oxygen dispersion and considers only the oxygen de- 
grees of freedom which couple directly to the Cu d x 2_ y 2 
orbitals. The results of the DCA calculation applied to 
the two-band Hamiltonian are presented in Section II I II 
The possible reduction of the two-band model to a single- 
band one, together with a detailed analysis of the single- 
band t-t'-U Hubbard model, is addressed in Section Hvl 
A discussion regarding the similarities and the differences 
between the two-band model and the single-band one is 
given in Section The conclusions of our study are 
reviewed in Section IVII 

II. FORMALISM 
A. The model Hamiltonian 

Band structure calculations^^! 2 !!, cluster calculation^, 
photoemissioni^ and other experiments show that the rel- 
evant Cu degrees of freedom are the d x 2_ y 2 orbitals which 
couple with the in-plane p x and p y O orbitals. All these 
degrees of freedom result in a five-band (four oxygen and 
one copper band) Hamiltonian in general. We have stud- 
ied the five-band model in detail 21 and have found that 
due to the strong Cu-0 hybridization, only the oxygen 
degrees of freedom which couple directly with Cu are rel- 
evant for the low-energy physics. Consequently, to a very 
good approximation, the five-band model can be reduced 
to a two-band one. 

The two-band model contains one Cu d x 2_ y 2 correlated 
band and one oxygen band. At every site the oxygen 
states are obtained by taking a linear combination with 
x 2 — y 2 symmetry of the four O p a orbitals which form 
a plaquette around the Cu ion. These are the only oxy- 
gen states which hybridize directly with Cu. However, 
it should be mentioned that these plaquette states are 
not orthogonal, two neighboring states sharing a com- 
mon oxygen atom. An orthogonal basis can be obtained 
by the procedure described in the original ZR paper 11 . 
First, applying a Fourier transform, translational invari- 
ant (Bloch) states are constructed. The Bloch states are 
orthogonal but not normalized, so they should be multi- 
plied by a normalization factor 

0(k) = [sin 2 (fc x /2) + sm 2 (k y /2}- 1 / 2 . (1) 

After normalization a complete and orthonormal set of 
oxygen states is obtained. 

In this basis the two-band Hubbard Hamiltonian can 
be written as: 

H = Y,E{k)cl (T c ka + E d d{ a d krT + V(k)(c{ a d ka + h.c) 
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FIG. 1: Coarse-graining of the Brillouin Zone in four cells 
(N c = 4) around K=(0, 0), (0, tt), (tt, 0) and (tt, tt). 

+ U n di in d ii (2) 

i 

We work in the hole representation and d\ a ( cj. CT ) repre- 
sents the creation operator of a Cu (O) hole with spin a 
and momentum k. The O band dispersion and the Cu-0 
hybridization are given by 

E(k) =E P - 8t pp f3 2 (k)sm 2 {k x /2)sm 2 {k y /2) (3) 

V(k) = 2t pd p~ 1 (k) (4) 

with t pp being the 0-0 hopping integral. The last term in 
Eq|5]represents the Coulomb repulsion between two holes 
on the same d orbital. We choose the commonly accepted 
values of the parameters, based on the band structure 
calculations of McMahan et al. — and Hybertsen et ali^ 
Because of the low density of oxygen holes (25% — 30%), 
we treat the Coulomb repulsion on O (given by U pp ) and 
the repulsion between nearest-neighbor Cu and O holes 
(given by U pd ) at the mean field level as a reasonable 
approximation. The effect will be an increases of our 
estimation for A = E p — E d by f7 p ^f + U pd (n d — n p ), 
where n d and n p are the average occupation of Cu and 
respectively O bands. A choice of U pp — 6 eV, U pd — 
1.3 eV and n p = 0.3 results in a increase of A by 1.3 eV. 

To conclude, we take in Eq. t pd = 1.3 eV, t pp = 
0.65 eV, A = 4.8 eV and U = 8.8 eV. 



B. DCA technique 

The DCA is an extension of the Dynamical Mean Field 
Theory (DMFT) 2 ^,. The DMFT maps the lattice problem 
to an impurity embedded self-consistently in a host and 
therefore neglects spatial correlations. The DCA maps 
the lattice to a finite-sized cluster embedded in a host. 
Non-local correlations up to the cluster size are treated 
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FIG. 2: Two-band Hubbard model DOS at 0% doping. 
The solid line is the d part of the DOS calculated at T = 
685 K. The value of the parameters is t P d = 1.3eV, t pp — 
0.65eV, A = 4.8eV, U = 8.8eV. The dashed line shows the 
DOS when t p d = 0. The chemical potential /i = 0. 



explicitly, while the physics on longer length scales is 
treated at the mean field level. This technique is ideal 
for the problems where short-ranged correlations are pre- 
dominant as in the high-T c materials. Here we calculate 
the properties of the embedded cluster with a Quantum 
Monte Carlo (QMC) algorithm. The cluster self-energy 
is used to calculate the properties of the host, and this 
procedure is repeated until a self-consistent convergent 
solution is reached. The self-energy and vertex functions 
of the cluster are then used to calculate lattice quanti- 
ties. A detailed description of the QMC-DCA algorithm 
is given in refi^ 

In the two-band model the oxygen degrees of freedom 
are not correlated, and therefore they are not included 
explicitly in the cluster. Their effect is fully contained 
in the cluster-host hybridization function and in the host 
Green's function. 

Here we consider a 2 x 2 cluster of Cu ions, which we 
believe to be large enough to capture the essential physics 
of Hubbard-type models. The 2x2 cluster will result in a 
coarse-graining of the BZ in four cells, as shown in Fig.^ 



III. TWO-BAND HUBBARD MODEL RESULTS 

The undoped materials have one hole per CuO unit. 
For t p d = the DOS is given by the dashed line in Fig. [21 
and the hole addition states would be of pure O char- 
acter. When the Cu-0 hybridization t pc t is switched on, 
the extra holes added to the oxygen band will scatter 
with the Cu spins and bound states will appear at the 
bottom of the oxygen band. This is illustrated by the 
solid line, which plots the partial d DOS which was ob- 
tained using the Maximum Entropy Method (MEM) 24 
for the analytic continuation of the QMC data to real 
frequencies. It can be noticed that the first hole addition 
states have a strongly mixed d and p character (the d 
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FIG. 3: Two-band Hubbard model phase diagram. 



character in the spectrum is large now) and an energy 
pushed well below the edge of the initial non-bonding 
oxygen band. Therefore only these states are relevant 
for the low-energy physics 3 ^. In the ZR theory these 
low-energy states which appear as a consequence of the 
strong Cu-0 hybridization are considered to be local sin- 
glets which move through the lattice like holes in an AFM 
background. Consequently the claim is that the physics 
can be described by a one-band t-J model. 

In order to determine the phase diagram we calculate 
a large number of susceptibilities which are relevant for 
spin, charge and superconducting ordering, both at the 
center and the corner of the BZ. For example, the Neel 
and SC critical temperatures, Tn and respectively T c in 
the phase diagram presented in Fig. |3| are determined 
from the divergence of the corresponding susceptibili- 
ties. The pseudogap crossover temperature T* is ob- 
tained from the maximum in the uniform magnetic sus- 
ceptibility when accompanied by a suppression of spec- 
tral weight in the DOS. Similar to what was found in the 
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SC for both electron and hole doped regimes. However 
the electron-hole symmetry is broken. In the electron 
doped case AFM persists to a much larger doping. On 
the contrary, SC disappears at a smaller critical doping 35 . 
These features of the phase diagram are in qualitative 
agreement with the experimental findings 2 . 

The one-electron spectral functions as measured with 
photoemission are also different. Our 2x2 cluster di- 
vides the BZ into four cells around K=(0, 0), (0, 7r), (tt, 0) 
and (tt, tt) (see Fig. ^| and approximates the lattice self- 
energy by a constant T,(K, lo) within a cell. Because of 
this coarse-graining, a comparison with ARPES is not 
possible apart from gross features. However, as the phase 
diagram shows, we believe that even our small cluster 
captures much of the physics of the cuprates. Here we 
want to stress the difference between the electron and 
the hole doped case within our 2x2 cluster approxi- 
mation. In Fig. 0] a) and b) we show the total d states 
DOS and the d coarse-grained K dependent DOS (which 
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FIG. 4: Total d DOS and coarse-grained K dependent d DOS 
at 5% doping, a) hole doping case, b) electron doping case. 



would correspond to the average over all k belonging to a 
coarse-graining cell of the single particle spectra A(k, uj)) 
for the hole and respectively for the electron doped case, 
at 5% doping. The total DOS looks qualitatively sim- 
ilar and at the chemical potential we see in both cases 
a depletion of states which indicates the presence of the 
pseudogap. The K dependent DOS is very different. The 
important feature which we want to stress is the location 
of the pseudogap in the BZ. In the hole doped case the 
pseudogap appears around (0, 7r). For the electron doped 
case we do not detect any suppression of states around 
(0, 7r) even though the pseudogap is clearly present in 
the total DOS. These features are in agreement with the 
photoemission experiments. The hole doped materials 
show Fermi pockets around (tt/2,tt/2) and gapped states 
around (0, ir)&. For the electron doped materials the pho- 
toemission spectral exhibit a gap near (ir/2,ir/2) and 
Fermi surface pockets around (0,7r). With the present 
cluster size the DCA cannot determine where in fc-space 
the pseudogap is, but it is interesting that it is not at 
(0, 7r). The presence of the pseudogap at (tt/2,tt/2) for 
the electron doped system can only be checked by in- 
creasing the cluster size, and this work is in progress. 

The electron and the hole doped susceptibility func- 
tions are also different both for the divergence temper- 
atures and the temperature and doping dependence. In 
Fig. [5] we show the uniform spin and charge suscepti- 
bilities versus temperature at 5% and 10% doping. A 
common feature for all cases is the existence of a char- 
acteristic temperature T* below which the spin response 
is suppressed and the charge one is enhanced. T* corre- 
sponds to the pseudogap (seen in the DOS) onset tem- 
perature. The suppression of the spin excitations be- 
low T* was also seen in NMRSi experiments and it was 
associated with the pseudogap. Besides these common 
features the electron and the hole doped susceptibilities 
behave differently. Generally, the maximum value of the 
spin susceptibility increases with hole filling. This means 



FIG. 5: Uniform Spill \spin (upper part) and charge \ch 
(lower part) susceptibilities versus temperature for different 
hole densities, n in the legend represents the number of holes 
per unit cell. 
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FIG. 6: Inverse of the d-wave pairing susceptibility Xsc versus 
temperature for different hole densities. 



that in the hole doped case, the spin susceptibility at the 
pseudogap temperature is strongly increasing with dop- 
ing unlike in the electron doped case where it decreases 
upon doping. At the same doping the hole doped spin 
susceptibility is much larger than the electron doped one. 
Another interesting feature is the very strong increase of 
the charge susceptibility for the electron doped case in 
the underdoped region (5% doping), suggesting a ten- 
dency towards phase separation^. 

Asymmetric behavior can also be noticed in Fig. [ft] 
where we plot the inverse of the d-wave paring suscepti- 
bility. Above T c the pairing susceptibility increases with 
doping in the electron doped case and remains more or 
less constant in the hole doped case. 

Because of the large Cu-0 hybridization the system is 
strongly covalent. For example in the undoped regime 



6 



0.8 


a) T = 538 k' ' 


0.7 






0.6 




• • nd 

■■■ n d o* n 




0.25 


_b) T = 538 K 

screened moment on Cu, T *Xi cal 


0.2 






0.15 






0.1 


■r 


i 



0.7 0.8 0.9 1 1.1 1.2 1.3 



n (filling) 

FIG. 7: a) The Cu occupation number, n<j, the unscreened 
Cu moment, fi 2 (Eq. versus hole filling, b) The screened 
Cu moment, T Xiocai (Eq. versus hole filling. 



where Si is Cu spin operator at site i, is shown. The 
main effect of the extra holes is to screen the spins on 
the Cu sites. The screening starts to be effective bel- 
low temperatures of about rj 0.5 eV (not shown). In 
the Zhang-Rice scenario an extra hole perfectly screens 
one spin on Cu forming a strongly bound on-site singlet 
which would contain a significant amount of the double 
occupied Cu configuration. So, our results do not contra- 
dict the ZR theory, but also do not exclude other scenar- 
ios where the extra holes form more complicated bound 
states which involve more than one Cu spin. Quantitative 
analysis based on the amount of screening as function of 
hole doping cannot give an answer to the validity of the 
ZR assumption, because, aside from the screening due to 
the oxygen holes, there are also non-local processes which 
contribute to the screening of Cu moments (for example, 
a possibility is the formation of inter-site spin singlets 
associated with the Resonance Valence Bond scenario). 



the Cu occupation number is only » 73%. The fact that 
the cuprates are strongly covalent was also observed in 
NMR measurements^. We notice that the system ex- 
hibits a slightly doping dependent covalency. This is 
shown in Fig.|7J-a where the Cu occupation number versus 
hole density is plotted. A constant covalency, equal to the 
one in the undoped regime (i.e. 0.73 Cu holes and 0.27 
O holes per site), would correspond to the dashed line. 
It can be noticed that, for the electron doped regime, the 
Cu hole occupation number is decreasing faster than the 
hole concentration, which indicates an increasing cova- 
lency with increasing electron doping. This happens be- 
cause at large electron doping, i.e. when the hole-filling 
of the CuOi plane is small, the effective hybridization 
is a result of a large V(k) in the BZ — . Increasing the 
number of holes, the BZ starts to fill up and a smaller 
V(k) will be responsible for the hybridization, and con- 
sequently the covalency decreases. For the hole doped 
regime, the extra holes go primarily on the oxygen band, 
and therefore we do not have a direct measure of the 
covalency. 

In Fig. [7\& the unscreened moment on the Cu orbitals 
is shown. It is defined as: 

M 2 = ((ridii - n dll ) 2 ) = n d - 2{n di i i n dii ) . (5) 

The difference between n d and /i 2 is a measure of the 
double occupancy with holes on Cu sites. In the electron 
doped regime the double occupancy is very small, but it 
increases substantially in the hole doped regime, which 
indicates that the low-energy hole addition states con- 
tain double occupied Cu configurations in a significant 
measure. 

In Fig. [3-b the screened moment on Cu, defined as 

T Xiocai = ^ ($~ ( T )^ + (°)> dT ' ( 6 ) 



IV. REDUCTION TO SINGLE-BAND 
HUBBARD MODEL 

Concluding that the electron-hole asymmetry is an in- 
trinsic property of the CuOi plane, we next address the 
cause of this asymmetry and the possible reduction to a 
one-band model. 

In Sec. IIIII we showed that, due to the Cu-0 hybridiza- 
tion, the addition of holes results in the formation of 
low-energy states, with an energy well bellow (« 1 eV) 
the initial oxygen band (see Fig. |2J). The reduction to 
a one-band model is based on the ZR claim that these 
states are singlets, i.e. spinless entities which can be 
regarded as holes moving in an antiferromagnetic back- 
ground. Due to the Monte Carlo nature of our calcu- 
lation, which does not provide a wave function for the 
ground state, we cannot determine directly the exact na- 
ture of these states. The most we can do is to compare 
the results of a two-band Hubbard model calculation with 
those of a one-band Hubbard model, and based on the 
similarities and differences that we might find to decide 
about the validity of the single-band approach. 



A. Zhang & Rice approximation and the derivation 
of the effective parameters 

In order to compare the two-band and the one-band 
models, we should first get an idea about the possible 
single-band model effective parameters. We discuss here 
two different approaches for calculating these parameters, 
both based on the assumption that the low-energy states 
are localized, and close to the ZR proposed singlets. 
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cell-perturbation 


U = 3.04 


J = 0.25 


J'«0 


t h = 0.477 


t e = -0.35 


£ J = 0.433 


t h ' = -0.03 


t e/ = -0.016 


t J ' = -0.003 


cluster calculation 




J = 0.192 


J' = 0.012 


t h = 0.452 


t e = -0.323 




t h ' -0.169 


t e ' = 0.078 





TABLE I: First row: parameters in eV calculated using cell-perturbation theory. Second row: parameters in eV calculated 
using cluster diagonalization. 



1. Cell-perturbation theory 

The cell-perturbation theory^ assumes that the ZR 
mapping is strictly true and therefore the low-energy 
states are genuine local singlets. Here and everywhere in 
the paper by local we refer to the oxygen orthogonal Wan- 
nier states which are different from the non-orthogonal 
plaquette states around the Cu ions. 

To deduce the one-band model parameters we work in 
the site representation. We can Fourier transform Eq. [21 
and write it as 

H = H Q + H 1 (7) 

where 

H = = Si ZA E ocld° + E d d\ a d ia 

i 

+V (cl a d la + h.c)) + Un dl ^n d ii ■ (8) 

Here i represents the site index. The oxygen operators 
Cj describe the orthogonal Wannier states. The ZR as- 
sumption implies that Hq is responsible for the forma- 
tion of the low-energy states (local singlets), and H\ will 
determine the hopping parameters. Therefore the cell- 
perturbation theory provides a means to determine the 
one-band parameters provided that the ZR assumption is 
correct. Elaborate calculations along this line were done 
inl& for a variety of multi-band parameters. In a first 
order approximation in Hi, the effective U is given by 

Ueff = E 2 +E° - 2E 1 (9) 

where E 2 , E 1 and E° represent the energies of the two 
(i.e. the ZR singlet), one (i.e. the bonding state) and 
respectively zero hole states of Eq.|SJ An important point 
is that Hi introduce three types of hoppings. If we denote 
with \2i >, 1 1, > and |0j > the lowest energy states of Hqi 
corresponding to two, one and respectively zero holes, we 
have the following hopping integrals 

t% = <2 l ,l ] \Hi\l l ,2 3 > (10) 
t% = < i ,l i |J?i|li,0 j > (11) 
t£ = < l^ljlHilOi^j > (12) 

t h (Eq. I1U[) describes the hopping of the ZR singlet, t e 
fEa. lllf) the hopping of the electron, and t J produces the 
exchange interaction 

J = At j2 /U eff . (13) 



The cell-perturbation theory applied to our model gives 
the parameters shown in the first row of Table. [I] 

We want to point out two things. First, the reduced 
Hamiltonian in the cell-perturbation theory is a t-t'-J 
model, 

H = -tY J b\b 3 -t> £ b\b 3 + J ]T SiSj , (14) 

<i,j> <<i,j>> <i,j> 

with different hopping parameters for the electron and 
the hole doped regions and with a value of the exchange 
interaction not determined by the quasiparticle's hopping 
(t h or t e ), but, as it is shown in Eq. ^]by t J . Therefore, 
a comparison with a one-band Hubbard model, should 
be done cautiously. Second, we want to stress that the 
value of the next-nearest-neighbor hopping terms (t' e and 
t' h ) is very small compared to the nearest-neighbor terms. 
The reason is that the initial oxygen-oxygen hybridiza- 
tion, t pp , results in an effective hopping term comparable 
in magnitude with the one resulting from the copper- 
oxygen hybridization, but with a different sign. This was 
also remarked ini&, and turns out to be an important 
observation for our final conclusions. 



2. Cluster calculation 

The other approach used for determining the parame- 
ters of the one-band model is based on a cluster calcu- 
lation. In order to estimate the nearest-neighbor hop- 
ping, the next-nearest-neighbor hopping and the ex- 
change terms, Eskes et a^considered two clusters, 
CuOy (which contains two nearest-neighbor Cu ions) 
and respectively CuOs (which contains two next-nearest- 
neighbor Cu ions). The exchange term is determined as 
the energy difference between the singlet and the triplet 
state of two holes on a cluster. For three holes on a clus- 
ter, the two energetically lowest states can be very well 
(98%) approximated with the bonding and anti-bonding 
states of a plaquette ZR singlet hopping between the two 
cells. Therefore, the differences between these two levels 
is two times the ZR singlet hopping, t h . In an analogous 
way, considering only one hole on a cluster, the electron 
hopping, t e , is determined. Using the cluster approach, 
our two-band model results in the effective parameters 
shown in the second row of the Table [I] 
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3. Comparison of the two approaches 

It can be immediately noticed that the two approaches 
produce different parameters, especially regarding the 
value of the next-nearest-neighbor hoppings. In the 
cluster calculation we obtain significant next-nearest- 
neighbor hoping terms, |t e '|/l ie | = °- 22 and \t h '\/\t h \ 
()..'! 7 with different sign for the hole and respectively elec- 
tron doped case (t h ' < 0, t el > 0). 

The reason for the discrepancy between the two ap- 
proaches is that, unlike the cell-perturbation method 
which considers local singlets, the cluster approach con- 
siders singlets between a Cu hole and a oxygen state 
formed on the plaquette around the Cu ion. Since the 
oxygen plaquette states are nonorthogonal it is possible 
to write them as a linear combination of many orthogonal 
oxygen states at different sites, i.e. the plaquette singlets 
are nonlocal states (in the orthogonal base) . It should be 
pointed out that at first glance this non-locality seems 
irrelevant (the overlap of the local oxygen states with 
the plaquette states—^ is about 96%), but apparently 
it turns out to influence the value of the next-nearest- 
neighbor hopping of the reduced Hamiltonian consider- 
ably. 

As pointed out by Emery and Reiter— , the plaquette 
singlets are in fact an admixture of ZR singlets and ZR 
triplets and this can result in a significant value of the 
spin-spin correlation on the oxygen sites^i. Of course, 
if the ZR singlet-ZR triplet admixture is significant, this 
will make a rigorous reduction to a one-band model im- 
possible, because now the low-energy states have also a 
spin component and therefore cannot be mapped into 
holes (spinless entities). 

It is worth pointing out that, in the cluster approach, 
the large value of the next-nearest-neighbor hopping 
terms results solely from the finite oxygen dispersion and 
the lack of hopping between the copper and the next- 
nearest neighbor oxygen plaquette state. On the other 
hand, in the cell-perturbation theory a copper - next- 
nearest neighbor oxygen hopping term is present. It re- 
sults in an effective next-nearest-neighbor hopping with 
a sign different than the one produced via oxygen-oxygen 
hopping. 



B. Possible reasons for the reduction to fail 

We believe that a comparison between the two-band 
Hubbard model and a single-band Hubbard model should 
be done with extreme caution. We want to stress the 
possible problems here. 

First, the reduction based on the ZR approximation, 
which results in a single-band t-J (or t-t'-J) model as- 
sumes the strong-coupling limit, i.e. a ratio U e ff/t ^> 8 
(the two-dimensional bandwidth is W — 8t). The low- 
energy density of states of the two-band model shown 
in Fig. |2 and Fig. 0] indicates a bandwidth of the or- 
der of the gap, showing that we are rather at the 
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FIG. 8: The relative double occupancy of the orbitals, 
< nfnj, > /n, versus hole filling, n, for different values of the 
ratio U jt of the single-band Hubbard model. 



intermediate-coupling than at strong-coupling. In the 
cell-perturbation theory we get U eff /t J = 7.02, which 
also suggests intermediate-coupling physics. Therefore, 
the question to be asked is whether the intermediate cou- 
pling regime, characterized by an effective repulsion of 
the same order of magnitude as the bandwidth, can still 
be well approximated by a second-order perturbation re- 
duced t-J model. 

Second, considering the previous objection, one may 
think that a reduction to the single-band Hubbard model 
in intermediate coupling regime, rather than to a t-J 
model, is more appropriate. However, serious problems 
arise from the fact that, in the ZR theory the nature of 
the antiferromagnetic correlations is different from that 
in the single-band Hubbard model, i.e. it is not directly 
related with the quasiparticle (ZR singlet or electron) 
hopping. Therefore, unless both the two-band Hubbard 
model and the one-band Hubbard model can be reduced 
to a t-J model, a comparison between them does not 
make much sense. Nevertheless, we believe that even 
when the effective repulsion is comparable with the band- 
width the second order perturbation theory which pro- 
duce the t-J model can be used successfully. We are 
going to discuss this at the beginning of the next section 
(Sec. HSU- 

Third, the non-locality of the low-energy states (in the 
sense discussed in Sec. lIV A~S|) can have very serious con- 
sequences beyond determining the value of the hopping 
parameters, making the single-band approach to fail com- 
pletely. 



C. t-t'-U Hubbard model results 

The t-J model results as a low-energy effective Hamil- 
tonian from the Hubbard model by projecting out the 
doubly occupied states. Therefore, the double occupancy 
of the site orbitals constitutes a measure of the validity 
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of this approximation. In Fig. [S] we plot the double occu- 
pancy of the site orbitals for different values of the ratio 
U/t. It can be noticed that for U/t > 6 the double occu- 
pancy is always smaller than 6%2&. This indicates that, 
even in the intermediate coupling regime the low-energy 
physics of the one-band Hubbard model can be well de- 
scribed by a t-J model. 

Even if, it is more natural to compare the two-band 
model with a t-t'-J (or a t-t'-J-J') model, this turns out 
from our perspective to be rather inconvenient due to the 
technical difficulties encountered by the QMC when ap- 
plied to such models. Therefore, we proceed by compar- 
ing the two-band model with a t-t'-U Hubbard model, 
focusing on the qualitative features rather than on a 
quantitative comparison. In the strong-coupling limit, 
the t-t'-U model reduces to a t-t'-J-J' model, with the 
constraint J' = J x (t' /t) 2 . Therefore, it is reasonable to 
assume that if the value of it! ft) 2 is not too large and the 
reduction of the two-band model to a single-band model 
is valid, the two models should exhibit similar physics. 

Assuming that the reduction to a one-band model in 
the spirit of the ZR approximation is possible, we should 
expect from Table. [I]the hopping parameters to be differ- 
ent in the hole doped and in the electron doped regions. 
On the other hand the exchange interaction, 



At 2 



(15) 



should be the same. 

Therefore, we study the single-band t-t'-U Hubbard 
model and address the following questions: 

1) How do the system properties depend on the ratio 
t/Jl 

2) What is the role of the next-nearest-neighbor hop- 
ing t'? 



1. The t/J dependence 

The values of the parameters in Table. U]show that in 
general the ratio \t/J\ is larger in the hole-doped regime 
than in the electron-doped case. In order to address 
the electron-hole asymmetry observed in the two-band 
model, in this section we study the properties of the 
single-band Hubbard model as a function of t/J, by keep- 
ing J (given by Eq.fHfl) constant and varying the hopping 
t. The next-nearest-neighbor hopping t' is set to zero. 

With respect to antiferromagnetism, with increasing 
t the Neel temperature at small doping and the critical 
doping where the antiferromagnetism disappear decrease. 
For example, at 5% doping, the antiferromagnetic suscep- 
tibility is diverging only for the small value of t shown in 
Fig. EH Assuming that the hole doped regime is charac- 
terized by a larger value of t/J, this feature is in agree- 
ment with the two-band model asymmetric behavior (see 
Fig.0. 
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FIG. 9: Antiferromagnetic susceptibility at 5% doping, for 
three different values of t, when J is constant. 
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The uniform spin susceptibility is shown in Fig. I1UI 
One can notice that an increase of t results in an in- 
crease of T* and a decrease of the spin susceptibility at 
T* . This together with the behavior of the susceptibility 
as a function of doping is in contrast to what was ob- 
served in the two-band model (see Fig.[3J) where the spin 
susceptibility is larger in the hole doped case and an in- 
crease (decrease) with doping of the susceptibility at T* 
for the hole (electron) doped regimes is found. 

The behavior of the d-wave pairing susceptibility as 
a function of t is shown in Fig. ^] The critical tem- 
perature increases with increasing t (the increase of T c is 
about 10% of the increase of t), as can be seen in the inset 
(a). This increase is much too large to be in agreement 
with the two-band model results even if, actually, for the 
two-band model we obtained a hole-doped T c larger with 
about 20 K than the electron-doped ona^S. By extrapo- 
lating the inverse of the ci-wave pairing susceptibility at 
28% doping (see inset (b)) it can be concluded that an 
increase of t results in an increase of the critical dop- 
ing where SC disappears. We also notice that, at small 
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FIG. 11: Inverse of d-wave pairing susceptibility versus tem- 
perature for different hole densities and hopping parameters. 
Inset -a) The critical temperature versus t at 5% (circle) and 
10% (squares) doping. Inset -b) Inverse of d-wave pairing sus- 
ceptibility versus temperature at 28% doping, for t — 0.37 eV 
(circles) and t = 0.52 eV (diamonds). 
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FIG. 12: Single-band t-U Hubbard model total and K depen- 
dent DOS at 5% doping. J = 0.22 eV, t = 0.45 eV 



FIG. 13: t-t'-U Hubbard model (solid line) and t-U Hubbard 
model (dashed line) phase diagrams for t = —0.45 eV, U = 
3.6 eV. For the t-t'-U Hubbard model t'/t = -0.3. 
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FIG. 14: t-t'-U Hubbard model. Uniform spin Xspin (upper 
part) and charge Xch (lower part) susceptibilities versus tem- 
perature for different hole densities. 



doping and above T c , a large t suppresses the pairing 
correlations. These features are in agreement with the 
asymmetry of the two-band model phase diagram. Nev- 
ertheless, we notice that, above T c and for both values of 
t, with increasing the doping the pairing correlations in- 
crease too. This behavior is characteristic in the electron 
doped regime of the two-band model, but cannot explain 
the hole doped regime where the pairing does not depend 
on the doping (see Fig. EJ) . The other difference between 
the two-band and the single-band Hubbard model is the 
value of the SC susceptibility critical exponent 7, which 
is much smaller in the two-band model case. 

The density of states and the K dependent DOS for 
the one-band t-U Hubbard model at 5% doping is shown 
in Fig. E| The one-particle spectra exhibit a pseudo- 
gap in the total DOS and in the K dependent DOS 
at (0, tt) point in BZ, similar to the hole doped spec- 
tra of the two-band Hubbard model. The single-band t- 
U Hubbard Hamiltonian is particle-hole symmetric and 



therefore cannot explain the one-particle spectra in the 
electron-doped regime of the two-band Hubbard model. 

At the end of this section we conclude the following: A 
single-band t-U Hubbard model (i.e. t' — 0) with a larger 
value of the hopping parameter for the hole doped regime 
cannot explain the electron-hole asymmetries observed in 
the two-band Hubbard model, especially the ones which 
characterize the one-particle spectral functions and the 
susceptibility functions. 

2. The t' dependence 

In this section we study the role of the next-nearest- 
neighbor hopping, t', in the single-band Hubbard model 

H = - t H h \ h i - l ' E b ^ + u E n w ■ ( 16 ) 

We choose the following parameters, U = 3.6 eV, t = 
—0.45 eV and t' = 0.15 eV. These parameters are close 
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FIG. 16: a) t — t' — U total DOS and coarse-grained K depen- 
dent DOS at 5% doping for t = -0.45 eV, t' /t = -0.3, U = 
3.6 eV. a) hole doping case, b) electron doping case. 



to the ones in Table. U resulting in J = 0.22 eV and 
J' = 0.02 eV. 

As for the two-band Hubbard model, we work in the 
hole representation, defined as the one where the filling 
1 + S corresponds to a hole doping S. Values of the fill- 
ing smaller than one correspond to the electron doped 
regime. We keep the sign of t' always positive. In or- 
der to avoid confusion we want to point out that in a 
t-J model the filling is always smaller than one. There- 
fore, in order to describe the electron and the hole doped 
regimes one has to employ the hole and, respectively, the 
electron representation. Accordingly, the sign of t' has to 
be chosen negative in the hole doped regime and positive 
in the electron doped case^t 

In Fig. EH the phase diagram of the t-t'-U model is 
shown with a solid line. For comparison, the phase di- 
agram of t-U Hubbard model (i.e. t' = case), which 
is symmetric with respect to hole and electron doping, 
is shown with dashed line. At half filling t' introduces 
an effective antiferromagnetic exchange J' = 4t' 2 /U be- 



tween the same sublattice spins and subsequently frus- 
trates the lattice. However at finite electron doping, t' 
favors the antiferromagnetism, making it persist up to a 
larger doping. On the other hand, in the hole doped case, 
the antiferromagnetism is always suppressed by t' . With 
respect to superconductivity, the presence of t' results in 
a smaller (larger) critical electron (hole) doping at which 
the superconductivity disappears. The asymmetry in- 
troduced by t' is in agreement with the one observed in 
the two-band model phase diagram. We find that t' has 
no major influence on the maximum superconductivity 
critical temperature T™ ax . 

The uniform spin and charge susceptibilities are shown 
in Fig.^J The spin susceptibility at the pseudogap tem- 
perature T* is strongly increasing with doping for the 
hole doped case and an opposite effect is seen for the elec- 
tron doped case. The downturn at T* in the spin suscep- 
tibility is much sharper for the hole doped regime indicat- 
ing a fast transition to the pseudogap physics. All these 
features are in very good qualitative agreement with the 
ones corresponding to the two-band Hubbard model. Be- 
cause of the similarity with the two-band model, it is also 
worth mentioning that in the electron-doped regime the 
charge susceptibility is strongly increased below T* in the 
undcrdoped region. 

The d-wave paring susceptibilities shown in Fig. II 51 ex- 
hibit asymmetric features, also in a qualitative agreement 
with those in the two-band model. In the electron doped 
regime, with increasing the doping, the pairing correla- 
tions above T c increase. In the hole-doped regime close 
to T c , the pairing correlations do not significantly depend 
on the doping. However, contrary to the two-band model 
behavior, at larger temperature an increase of pairing 
correlations with doping is observed. The magnitude of 
this increase is smaller than in the electron-doped case 
and a larger value of t' (e.g. t' « 0.4£, not shown) will 
reduce it further, improving the resemblance with the 
two-band model. 

In Fig. UHlwe present the DOS of the t-t'-U Hubbard 
model at 5% doping. The one-particle spectral functions 
resemble the corresponding two-band Hubbard model 
ones. The presence of the t' parameter is responsible 
for the location of the pseudogap in the BZ. 

The necessity of the t' in explaining the measured 
ARPES line shape and the electron-hole asymmetry was 
realized early on^iS. Representing hoppings in the same 
sublattice, this parameter is not severely renormalizcd 
by the AFM background, and consequently its influence 
turns out to be important. Exact diagonalization re- 
sults^ of a t-t'-J model are in agreement with ours. The 
t' hopping process lowers the kinetic energy and moves 
the quasiparticle position from (tt/2, ir/2) to (0, tt) in the 
electron doped case . The Neel-like configurations, which 
do not hinder this process, are stabilized. In the hole 
doped case the t' hopping does not lower the kinetic en- 
ergy of quasiparticles and it is not energetically favorable, 
therefore leading to a suppression of AFM at all dopings. 

The main conclusion of this section is that a one-band 
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t-t'-U Hubbard model describes qualitatively well the 
physics (i.e. the phase diagram, the one-particle spectra 
and the two-particle response functions) of the two-band 
Hubbard model, provided a significant value of the next- 
nearest neighbor hopping (i'/t rs 0.2 — 0.5) is considered. 
However, besides all these similarities there are also some 
important differences which we are going to emphasize in 
the next section. 



V. DISCUSSION 

In general, the deduction of an effective low-energy 
Hamiltonian implies two steps. First, defining the low- 
energy states and second, projecting the resolvent op- 
erator, G(E) — (E — -ff) -1 , on the subspace spanned 
by these low-energy states^. The inverse of the pro- 
jected operator can be written as E — H e ff(E), where 
H e ff is the low-energy Hamiltonianii. This procedure 
is equivalent to finding an Hamiltonian which produces 
the same one-particle, two-particle, three-particle, etc., 
spectral functions on the energy range considered to be 
"low-energy" . 

Rigorously, in order to prove that the one-band model 
is the effective Hamiltonian which describes the two-band 
Hubbard model low-energy physics we should compare 
not only the one-particle and the two-particle spectra, 
but also all higher order correlation functions. How- 
ever, we believe that the comparison of only the one- 
particle and the two-particle spectral functions is com- 
pelling enough, especially since the experimental infor- 
mation is also obtained by measuring the response func- 
tions behavior (and in almost all cases the two-particle 
operators or the one-particle ones, as in photoemission, 
are involved). It is also true that a comparison of the 
dynamic susceptibilities would be required, but with our 
Quantum Monte Carlo based algorithm the calculation 
of these quantities for the two-band model is extremely 
computational resources consuming and has not been 
done yet. However partial information about the relevant 
excited states is contained in the temperature behavior 
of the static susceptibilities. 

The main conclusion of Sec. IIVI is that a t-t'-U Hub- 
bard model describes qualitatively well the physics of 
the two-band Hubbard model, but only if a substantial 
next-nearest-neighbor hopping is considered. However, 
the calculation in Sec. IIV A II (first row of Table. P) and 
the more rigorous results by Feiner et alm&, show that 
in a strict ZR picture the next-nearest-neighbor hopping 
is negligible. Therefore it is difficult to explain the two- 
band Hubbard model physics assuming the formation of 
local ZR singlets. For hole doped systems, a significant 
value of t' can be obtained only if the extra holes form 
nonlocal bound states with the existing Cu holes, pre- 
sumably something close to the plaquette singlets. Of 
course we have no reasons to discard other states spread 
over even more oxygen sites, which can result in a mag- 
nitude of the hopping parameters different (probably not 



too much) form the one obtained by cluster calculation 
(second row of Tabic. Q). In the electron doped systems, 
the doping dependent covalency shown in Fig.|7J-a, clearly 
indicates that the hybridization of the Cu with the O 
states at different sites is important. A doping dependent 
covalency should also imply doping dependent parame- 
ters. 

The cluster calculation which allows the formation of 
non-local (plaquette) low-energy states, unlike the cell- 
perturbation (or strict ZR) approach, provides a value 
of the hopping parameters which captures qualitatively 
the physics of the two-band model. However, we do 
not believe that finding the exact value of the one-band 
Hubbard model parameters is a relevant or even a well 
addressed problem, because the non-locality of the low- 
energy states implies that the two models are not equiv- 
alent, as was pointed out by Emery et alw^. Aside from 
the similarities between the two-band and the t-t'-U Hub- 
bard models discussed in Sec IIV C~2l we also find some 
differences. 

For example, one important difference can be observed 
in the d-wave pairing susceptibility (Fig. I15| and Fig.[SJ). 
In the two-band Hubbard model the critical exponent 7 
which defines the divergence of the susceptibility at T c , 
is much smaller (around w 0.4 at finite hole doping) than 
the one characteristic to the one-band model (around » 
0.6), indicating larger fluctuations^. 

Both the cell-perturbation and cluster calculation pro- 
vide a larger nearest-neighbor hopping t for the hole 
doped region. According to the analysis presented in 
Sec. IIV C II this should result in both larger T* and T c . 
However the two-band model results does not indicate 
that this is the case, the respective critical temperatures 
being not very different in the electron and hole doped 
regimes. 

Based on our comparison we can draw the following 
conclusions. The one-band Hubbard model retains much 
of the two-band Hubbard model physics, but a signifi- 
cant next-nearest-neighbor hopping (t'/t w 0.3) should 
be provided. If the purpose of the investigation is the 
study of the basic physics like the SC mechanism, the 
proximity of AFM, SC and pseudogap we believe that a 
one-band t-t'-U Hubbard model should be good enough. 
On the other hand, if the purpose is to describe more sub- 
tle features like the ones which may result from the finite 
value of the spin correlation on oxygen, or if a quanti- 
tative material specific calculation is desired, the single- 
band model approach fails. Obviously also the single- 
band model should not be used to describe spectral fea- 
tures at energies above 0.5 eV such as the optical, elec- 
tron energy loss and inelastic x-ray scattering results. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we use the DCA to calculate the prop- 
erties of the two-band Hubbard model. The 2x2-site 
cluster phase diagram resembles the generic phase dia- 
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gram of the cuprates and exhibits electron-hole asymme- 
try. We also find asymmetric features for the one-particle 
spectral functions and for the relevant susceptibility func- 
tions. These characteristics are in qualitative agreement 
with experimental findings. 

We address the validity of the single-band Hamiltonian 
as the effective low-energy model for the cuprates. We 
discuss the possible problems which may cause the fail- 
ure of the reduction from two-band to one-band and also 
show that, depending on the approximations involved, 
the value of the one-band Hubbard parameters (espe- 
cially the next-nearest-neighbor hopping) can be signifi- 
cantly different. 

We use DCA to study the role of the different param- 
eters in the single-band t-t'-U Hubbard model and com- 
pare the phase diagram, the one-particle and the two- 
particle response functions with the corresponding two- 
band Hubbard model ones. We conclude that the two 
models exhibit similar low-energy physics provided that 
a significant next-nearest-neighbor hopping t' is consid- 
ered. The parameter t' is also the main culprit for the 
electron-hole asymmetry of the cuprates. 

The large value of t' needed for a qualitative agreement 
between the two models cannot be obtained in a strict 
ZR picture where the extra holes form local singlets with 
the existing Cu holes. Plaquette singlets, which in the 
oxygen Wannier representation are not local, and pre- 
sumably other spatially extended states can provide a 
larger value of t'. The doping dependent covalency in 
the electron doped case also indicates that the non-local 
Cu-0 hybridization is important. However, the forma- 
tion of non-local low-energy states also implies that they 
are not real singlets and consequently cannot be rigor- 
ously mapped into holes and therefore the two models 
are not equivalent. 

We also point out some differences between the two 
models. In the two-band Hubbard model the fluctua- 
tions in the d-wave pairing channel above T c is much 
stronger. The deduction of the parameters both in cell- 
perturbation and cluster approach results in a larger 



nearest-neighbor hopping t for the hole doped regime. 
However the critical temperatures T* and T c in the two- 
band Hubbard model are approximatively the same in 
both regimes, quite different from what should be ex- 
pected. 

The conclusion is that a single-band Hubbard model 
with a significant value of the next-nearest-neighbor hop- 
ping (t'/t w 0.3), captures the basic physics of the two- 
band Hubbard model, including the proximity of antifer- 
romagnetism, superconductivity and pseudogap and ex- 
plaining the electron-hole asymmetry seen in the phase 
diagram, one-particle and two-particle spectral functions. 
However, the single-band Hubbard model is not entirely 
equivalent to the two-band Hubbard one and we believe 
that it is not suitable for quantitative material specific 
studies or for describing more subtle features which may 
results from the non-locality of the low-energy states. It 
is also not suitable to describe physics which implies ex- 
citations with energy scales larger than w 0.5 eV. 
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